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Background Of The Invention 

This invention relates to methods and apparatus for pricing financial derivatives. 

Current pricing models for financial derivatives are based on two pricing methods: (1) 
the "Classic Arbitrage Argument 5 ' (CAA), or (2) "Monte-Carlo" method. The CAA states the 
return on a perfectly-hedged investment portfolio is precisely equal to the return on U.S. 
Treasury bonds. Under the CAA method, it is assumed that investors are capable of 
continuously hedging their portfolios in infinitesimal increments with no costs. Moreover, the 
CAA method assumes that the parameters of the model do not change during the lifetime of the 
portfolio and that the markets will always precisely follow the pricing model. 

Pricing models based on the CAA method are often referred to as "Risk Neutral Pricing" 
(RNP) models, because no risks are assumed by the CAA method. Under the CAA method, 
management of all risks and all discrepancies between the pricing model and the real market is 
passed to the human users of the model. RNP models are implemented for production in a form 
of numerical solutions to partial differential equations via finite difference schemes, finite 
element methods, closed form formulas, or by various tree type methods, including binomial and 
trinomial ones. For example, the Black-Scholes pricing method is an example of a method that 
is based on the CAA method and idealistic RNP assumptions. 

Monte Carlo pricing models, on the other hand, are based on statistical and probabilistic 
evaluation of possible future price scenarios. These scenarios are typically generated with the 
aid of random number generators. Evaluation of derivative prices is modeled via the CAA, 
discounted payoff method, or a combination of both. In either case, the value of the risk does not 
explicitly enter the evaluation formula. 

Therefore, RNP pricing models do not include market risks and costs in the valuation of a 
financial derivative. As a result, these pricing models generate valuations of financial 
derivatives that are not completely consistent with actual market behavior. Although the risks 
have been reflected in the market prices by various non-mathematical methods, the lack of 



appropriate mathematical modeling has precluded the explicit inclusion of risks into these 
pricing and hedging models. 
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Summary o f Invention 

The present invention overcomes the limitation of the prior art by introducing the risk 
premiums into a pricing model for calculating values of derivative securities. 

It is an object of the present invention to provide a pricing model for calculating values of 
derivative securities that consider various risk premiums. 

It is a further object of the present invention to provide a pricing model for calculating 
values of derivative securities that is consistent with market behavior. 

One or more of the foregoing objects are met in whole or in part by the present risk- 
adjusted pricing method and computer apparatus. The method includes determining the 
underlying security and derivative, determining the risks and trading costs associated with the 
derivative, formulating the RAP equation for the derivative, solving the RAP equation for the 
derivative and outputting a value of the derivative based on the solution of the RAP equation. 
The computer apparatus employs this method in a financial analysis application. The system 
includes an input unit for taking inputs of derivative characteristics and model parameters, a 
processing unit for computing a value of the derivative based on the input characteristics and 
parameters using a risk-adjusted pricing equation, and an output unit for displaying the value of 
the derivative. 



Description of the Drawings 
Figure 1 is a flowchart of a preferred embodiment of the RAP evaluation software. 
Figure 2 is an output generated by use of a preferred embodiment of the RAP evaluation 
5 software that describes the Black-Scholes implied volatility. 

Figure 3 is an output generated by use of a preferred embodiment of the RAP evaluation 
software that describes the difference between the Black-Scholes call price and the risk-adjusted 
call price. 

Figure 4 is an output generated by use of a preferred embodiment of the RAP evaluation 
;dlO software that describes the difference between the Black-Scholes delta and the risk-adjusted 
;if delta. 

Jf 5 ? Figure 5 is an output generated by use of a preferred embodiment of the RAP evaluation 

software that describes the difference between the Black-Scholes gamma and the risk-adjusted 
«_ gamma. 

p5 Figure 6 is the source code of a preferred embodiment of the RAP evaluation software 

M written in C/C++ language. 

,jg Figure 7 the article by Milan Kratka, Risk, April 1998, referenced in the specification. 
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Description of the Preferred Embodiment 

The present invention overcomes the limitation of the prior art by introducing the risk 
premiums into a pricing model for calculating values of derivative securities. A derivative 
security is a security that derives its value from another underlying security, commodity, or 
object. Derivative securities include futures, forwards, options on stock, options on index, 
options on futures, basket options, swaps, swaptions, and others. 

Contrary to the current RNP models based on the CAA or Monte-Carlo, the present 
invention is based on the Advanced Arbitrage Argument (AAA), which incorporates the value of 
the various risks associated with derivatives into the valuation of those derivatives. The AAA 
states that the expected return on an investment should not only equal the return on the U.S. 
government bonds (as assumed by CAA), but should also include expected investment costs 
(such as ask bid spreads, brokerage fees, borrowing fees, etc.) and premiums for taking 
investment risks (such as implied volatility risk, positive portfolio variance risk, price gap risk, 
etc.). The AAA can be summarized by the following template pricing equation, which is valid at 
any step of time during the life of an investment: 
(1) 

Expected Return Return Rate of Expected Rate of Expected 

on = on + Investment 4- Risk 

Investment Government Bonds Costs Premium 

Return on government bonds is often referred to as the "risk-free rate." Expected risk 
premium compensates the investor for taking an exposure to various risks associated with the 
investment, similar in function and form to insurance premiums. The present invention 
implements the AAA through Risk Adjusted Pricing (RAP) methodology, which provides a 
means to quantitatively enumerate individual risk premiums into the pricing equation. 

Because the AAA is valid at any step of time during the life of an investment, the RAP 
methodology can formulate the expected return on investment, rate of expected investment costs, 
and the rate of expected risk premiums at any given time in the future as a function of price, 



time, sensitivity, and other parameters. 

Suppose that an investor prices an option O on an underlying security S, which follows 
the typical lognormal process with the drift (or the trend) ft and the volatility a : 

(2) dS = juSdt + aSdz 

5 Here the term dt represents a small increment of time, dz represents a normalized stochastic 
process, and dS represents a change of the price of the underlying security over this time interval 
du By forming a portfolio P consisting of a unit of the option O and - 8 units of the underlying 
security S, the portfolio P can be expressed as 

(3) P = 0-SS 

HiO and is subject to the AAA. The amount 8 (delta) is usually called the hedge. 
J In the current RNP models, this portfolio's delta is adjusted continuously and without 

M costs. In the present invention, however, it may be adjusted at discrete steps of time based on a 
B re-hedging strategy of portfolio P. Though transaction costs occur at these discrete steps of time, 
the present invention spreads the costs over the entire time interval between re-hedging points, 
|15 similar to the customary practice with the dividend rate, thus defining a rate of transaction costs, 
,m rrc- These costs may include, but are not limited to, ask bid spread, slippage, brokerage and 
q exchange fees, short sales taxes, etc. There may be other trading costs associated with the 
maintenance of the portfolio. For example, short-selling of £ units of stock accrues fees for 
borrowing the stock in a rate r B , resulting in a cost r B 8Sdt over a small time interval dt. 
20 Similarly, a short position in a dividend paying stock or an index with dividend yield r D results in 
the costs r D 8 Sdt over a small time interval dt. 

The invested amount in this portfolio could have been invested in the government bonds 
resulting in a riskless return r. For portfolio P consisting of 8 units of stock S and option O on 
this stock, this amount may be equal to 0-8S . For portfolio P consisting of 8 units of futures S 
25 and an option O on this futures, this amount may be equal to just 0. Because the portfolio P 
cannot be hedged continuously, its value has positive variance, and thus cannot be considered 
riskless. It is subject to the positive variance risk: the higher the portfolio variance the higher the 



risk. The investment should thus generate a risk premium tvar proportional to this variance risk. 

Besides the positive variance risk, the investment is also exposed to risks that the market 
will not follow our assumptions. A possibility of an option price O jump without a change in the 
underlying value S generates implied volatility risk. The investment should thus generate a risk 
5 premium rate r a proportional to this implied volatility risk. A possibility of a price gap in the 

underlying security S generates the gap risk. The investment should therefore generate a risk 
premium rate yqap proportional to this gap risk. In the presence of other considered risks, for 
instance: credit risks, interest rate risks, liquidity risks, daily price limit risks, etc., we would 
include the corresponding risk premiums. The value of the risk premium is not restricted to the 

30 positive values. A negative risk premium corresponds to a negative risk exposure. 

u In the present invention, the expected return E[dP] on the investment in the portfolio P 

jlJ over a short time interval dt can be expressed as: 

¥t (4) E[dP] = E[dO-SdS] 

Following the AAA the "Risk-Adjusted Pricing" (RAP) template equation for an option on a 

i§5 stock or an index S is: 

3 ^ E[dO - SdS] = [r(0 - SS)dt]+ [r TC dt + r B SSdt + r D SSdt]+[r VAR dt + r a dt + r GAP dt\ 

% where the rate of borrowing fees r B and the dividend yield r D are zero for negative values of S . 
Similarly, the RAP template equation for an option on a futures contract S can be written as: 

(6) E[dO - SdS] = [rOdt]+ [r TC dt]+ [r VAR dt + rjt + r GAP dt\ 

20 From the assumption about the underlying security price process: dS = juSdt + oSdz and from the 
Ito's lemma of stochastic calculus, the stochastic differential equation is: 

(7) dO = d t Odt + \cr 2 S 2 d 2 ss Odt + d s OdS 

Re-hedging of the portfolio P over a time interval dt to the value S = d s O makes the expected 
value of E[dO - SdS] equal to: 
25 (8) E[d t Odt + ±a 2 S 2 d\ s Odt + d s OdS-Ms]=d t Odt+ 

and at the same time it minimizes the variance of the portfolio. Re-hedging consistently to a 
different delta S would make the expectation term E[dO - SdS] equal to: 
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(9) d t Odt + ± a 2 S 2 d 2 ss Odt + (d s O - S)fiSdt , 

and it would increase the portfolio variance by (d s O-S) 2 cr 2 S 2 dt . For investors with strong 

trend expectations, it is thus optimal to re-hedge to a trend-adjusted delta that differs from the 
value d s O. 

5 In general, it is difficult to estimate the future trend term ju from the market. It is a 

subjective issue. If the trend was known in advance and had a significant value, the underlying 
price would have been already adjusted accordingly. There is much uncertainty about trends. 
Each time somebody buys a security, there is somebody who sells it. Similarly, some investors 
might under-hedge, while others might over-hedge. For the pricing purposes we could assume 
□0 that the portfolio P is being hedged with the delta given by the value d s O. 

1J The rate of the transaction costs r TC related to the re-hedging of the portfolio typically 

aJ depend on the momentary market conditions that include the size of the ask bid spread, liquidity 
B of the market, and size of the trade. Since large orders tend to move the market, typically the 
transaction costs per unit contract will be a non-increasing function of the total size of the trade. 
Sf 5 If the function C denotes an aggregate trading cost per contract as a function of the trading size 
II for the individual market, then the rate of the transaction costs between two re-hedging points 

% over the time interval A, can be estimated as: 

c(n\AS\)aS 

(10) r TC = v 1 " , 

where the expected actual size E[AS] of the re-hedging trade, i.e. the number of traded units of 
20 the security S, is given by the formula: 

that can be found in the article by Milan Kratka, Risk, April 1998, at 67-71, shown in Figure 7, 

which is herein incorporated by reference in its entirety. Assuming that the prevailing rate of 

transaction costs in the market is constant, in the present invention the rate of the transaction cost 

25 is equal to: 

(12) r K = C^aS\T\-^. 
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Notably, the longer the time between re-hedging A t , the lower the trading costs rate rrc- 

To estimate the rate of the risk premium r VAR compensating the investor for taking the 
positive variance risk, note that the higher the variance the higher the risk premium. A variance 
of a normal process increases linearly in time, and the rate of variance risk premium can be 
5 modeled by a linear function of the portfolio variance: 
(13) r VAR A t =R*vn(M>\ 

where R is the variance risk coefficient, which can also be viewed as a risk aversion coefficient. 
The higher it is, the more risk premium should be collected for the variance risk. The variance of 
the portfolio var(AP) between two re-hedging points with rehedging to d s O over an interval A, 

4 0 can be mathematically estimated as: 

tJ (14) var(AP) = 3a 4 S 4 T 2 A 2 t , 

ij where T = d 2 ss O as it is standard in the industry. Since the adjustment to the portfolio is done in 
M discrete steps, measure of the variance of the portfolio between adjustment may be subject to a 
definition discussion. Nevertheless, the dimensional analysis of the expected formula for the 
35 variance confirms the proportionality of the variance to <x 4 S 4 r 2 A 2 , . Therefore in the present 

I invention the rate of the variance risk premium is: 

I (15) r VAR = 3R**S 4 r 2 A r 

Notably, the longer the time between re-hedging A t , the higher the variance, and the higher the 
rate of the variance risk premium r VA R. Since the variance of the hedged portfolio P should be 

20 less than a variance of its individual components, an upper bound for the variance risk premium 
from the variance of the underlying security is: 
(16) r VAR <RS 2 a 2 S 2 . 

This estimate is actual typically close to expiration of the option 0, when T becomes large. 

The risk premium related to the implied volatility risk exposure depends on a probability 
25 that the implied volatility a may change, and it depends on the sensitivity of the portfolio value 
to this volatility change. The sensitivity of the portfolio value P to the change in a is equal the 
sensitivity of the option value O to the change in a . Thus, in the present invention, the implied 
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volatility risk premium rate as: 

(17) r CT =v3 CT 6>, 

where v is the coefficient of dependency. 

The risk premium tqap for the underlying price gap exposure should be proportional to 
5 the probability of the price gap, as well as to the size of the price gap. Due to the Ito's lemma, 
the change of the portfolio value due to the price gap AS is equal to: 

(18) ap = \t(asJ, 

the gap risk premium can be modeled as: 

(19) r GAP =g\a 2 S 2 T, 
is d?0 where the expected price gap size is scaled. 

|^ If the financial markets are efficient, then the market pricing will reflect the market value 

W of risk, average prevailing trading costs, as well as competitive re-hedging strategies. There is 
ffl never a guarantee that the free market will value options in any particular way. But we can 
assume that it does, build a pricing model and use it to verify our assumptions. The optimal re- 
J35 hedging strategy in the market dictates minimization of the transaction costs rjc and the premium 
-u for the positive portfolio variance war- This leads to a simple algebraic optimization task 
,S yielding the solution in terms of optimal re-hedging time interval: 

(20) A,=A, 

<r 2 S 2 P 

where the coefficient This optimal re-hedging time interval corresponds to an 

20 optimal re-hedging strategy defined by re-hedging with every change of: 

_\_ 

(21) k s =±oSjK t =±k-T l 
in the underlying security S. 

This optimal strategy balances prevailing trading costs rjc against the positive variance 
risk premium rate r V AR- Taking into account that the variance estimate is bounded from above by 
25 the variance of the underlying security component of the portfolio, in the present invention the 
quantitative enumeration of all individual components in the template RAP equation for the 
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option on a stock or an index can be summarized as: 

(22) 0 + ^a 2 S 2 TM = r(0 - SS)+ (r B +r D )5S + vV, 

where the RAP volatility adjustment coefficient is denoted by: 

( i "\ 



(23) M = l + g-min 



mT\2RS 2 r 



2t-i 

J 



dO c)^) 
5 the time decay sensitivity theta is denoted by 0 = — , the hedging sensitivity delta by S = — , 

dt dS 

d 2 O dO 
the gamma sensitivity r = — - , and volatility sensitivity by vega by V = — . The coefficient M 

dS da 

shall always be positive to preserve the phenomenological interpretation of the RAP equation. 

Note that the borrowing fees rate r B is zero when delta is negative. Similarly, in the present 
f. invention, the RAP equation for the option on a futures contracts can be written as: 
% (24) © + ^a 2 S 2 rM = rO + vV. 

]i These RAP equations are partial differential equations that can be solved by using 

various available numerical methods. These RAP equations are typically coupled by so called 
J boundary conditions. These boundary conditions usually arise as a result of the option 
~f specifications. Similarly, numerical methods are used for handling extra conditions that arise 
^5 from the options specifications. For instance when dealing with an early exercise, the numerical 
solution to the RAP equation will make sure that the value of the option at any particular time 
does not fall below its exercise value during typical market conditions. 

The model coefficients in the RAP equations, such as r, a , or m, do not need to be 
constant, but for instance they may be functions of the time t. The RAP methodology includes as 
20 a special case the idealistic RNP models. For example, if the possibility of continuous costless 
hedging without presence of any risks is assumed, the AAA becomes same as the CAA, and the 
RAP model becomes identical to the Black-Scholes model. On the other hand, the 
phenomenological explanation and tangible modeling of individual risks and costs included in 
the RAP methodology allows the RAP models to better explain market pricing even in changing 
25 or volatile market environment. The RAP methodology offers not only more precise derivatives 
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pricing tools, but also offers more advanced risk management tools through its computation of 
the derivatives sensitivities. The RAP methodology offers a direct way to extract more specific 
information about value of the market risks and prevailing costs via the inverse problem of 
fitting the RAP model parameters with sufficient amount of the market data, which can be very 
5 useful for the better risk monitoring and portfolio management. 

Figure 1 shows a flow chart of the preferred embodiment of RAP evaluation software that 
computes the value of the derivative by solving the respective RAP equation for that particular 
derivative. At block 10, derivative characteristics are input for the particular derivative whose 
value is being calculated. For example, if the derivative is an option, the derivative 

d;0 characteristics may include the type of the option (vanilla puts or calls, spreads, or exotics), the 
payoff profile, time to expiration, and early exercise feature (American, where the right to 
exercise the option applies at any point in time until maturity of the option, European, where the 

B right to exercise the option applies only at maturity of the option, or Bermudian, where the right 
to exercise the option applies only at specific points in time). These characteristics are input 

35 manually by the user or automatically by an application software, such as Infinity™ , 

1 PCQuote™, or MicroHedge™. 

S At block 15, model parameters, like derivative characteristics, are input manually by the 

user or automatically by application software. In contrast to derivative characteristics, which 
describe the particular derivative, model parameters describe the behavior of that particular 

20 derivative in given markets (equity market, interest rate market, commodity market, foreign 
exchange market, etc.). For an option, the model parameters could include volatility, risk 
premium coefficients, trading costs, and risk-free rate. 

At block 20, the computation precision is established by setting the resolution in terms of 
the number of time-steps from the maturity time to the current time, and the number of price- 

25 steps and the lower and upper bounds for the underlying price. The lower the resolution, the 
faster the computation. The higher the resolution, the more precise the computation. These 
functions described in blocks 10, 15, and 20 can be implemented as an input unit. 
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Once the derivative characteristics, model parameters, and computation precision are 
input, the value of the derivative at maturity is initialized by the payoff function of this derivative 
at block 25. For example, the payoff function of the plain vanilla call option is equal to the 
maximum value of zero and the difference between the underlying security price at the maturity 
5 time and the strike price of the call option. By initializing the value of the derivative at maturity, 
this necessarily establishes the maturity time at block 30, which is subsequently used in block 35 
with the input time-step from block 20 to establish the new time at which the value of the 
derivative will be computed using the RAP equation for that specific derivative. In block 40, the 
new time is compared to the current time, and if the new time is greater than or equal to the 

do current time, then the value of the derivative is calculated in block 45 at the new time using the 
RAP equation for that respective derivative. 

^ After the value of the derivative is calculated, the software returns to block 35 to 

B calculate a new time. That new time is then again compared to the current time in block 40 as 
explained above. When the new time is less than the current time, the value of the derivative is 

J35 computed up to the current time and displayed as an output at 50. The functions described in 

I blocks 25, 30, 35, 40, and 45 can be implemented as a processing module and block 50 can be 

n implemented as an output unit. 

In block 45, in the preferred embodiment, the RAP equation is solved by the semi- 
implicit Crank-Nicolson finite difference scheme with the explicit method used for the nonlinear 

20 terms in the RAP equation. The RAP equation may be solved also by using other finite 
difference schemes, finite element, Runge-Kutta ordinary differential equations, or other 
numerical methods. In block 45, not only are the derivative's values computed, but also its 
partial derivatives. 

Figure 6 shows the source code for a preferred embodiment of the RAP evaluation 
25 software, herein incorporated by reference. The software is useful for evaluating and pricing 
financial derivatives with use of a computer or any other electronic calculation aid and allows 
entry by a user of various inputs, while producing a tangible price output. Some preferred 
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outputs of such a device include a discrete price or value of a derivative, graphical 
representations of price, as shown in Figures 2-5, and reports of values of derivatives. Of course, 
the output can be modified, added to, or supplemented in myriad ways depending on the user's 
preferences. Further, as stated above, the method and device of the present invention allows the 
user to input discrete risk factors and utilize those risk factors in any ways preferred by the user 
to produce the above-mentioned results. 

The graphical representations shown in Figures 2-5 can be generated by using the values 
calculated for the derivative by the solution of the RAP equation and entering the values into a 
commercial graphing package, such as Microsoft™ EXCEL™. 

With reference to Figure 6, an exemplary computer code listing has been provided 
commensurate with the disclosed invention. The source code listing of Figure 6 is written in 
C/C++ and is best utilized with respect to a call or put security based on a European variety 
option. The code computes the value of a security based on a time reference as shown in the 
flow chart of Figure 1, above. The computer code can be run on any well known computer 
system as is well known to those of ordinary skill in the art. The code first requests the 
following inputs: 1) T, time to maturity in years; 2) K, strike price in dollars; 3) S, current price 
of the underlying security in dollars; 4) Y, option type (put or call); 5) sig, implied underlying 
volatility; 6) mu, expected trend rate; 7), R, variance risk premium; 8) C, trading cost/slippage 
coefficient; 9) v, volatility risk premium coefficient; 10) g, gap risk coefficient; 11) r, risk free 
rate; 12) rb, borrowing fees rate; 13) rd, dividend yield rate; 14) nt, number of time steps; 15) ns, 
number of price steps; 16) Is, lower bound of underlying price; and 16) us, upper bound of 
underlying price. The program initializes the variables, and creates the appropriate array of 
matrix coefficients. The program then verifies that the coefficients are within their normal 
ranges. The program then processes the variables and numerical inputs using the RAP method, 
as stated in the comments in the program itself, as shown. The program then stores the two- 
dimensional array as computed, and produces a value of the security. The program steps also 
compute and produce required theoretical values of securities in real time, and display or report 
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on the same. Further, use of the computer program model and steps allow real-time reporting 
and monitoring of the value of the underlying risks, to extract the value of the risk embedded in 
the security, to generate and display a delta hedging parameter (sensitivity of the option value to 
change in the value of the underlying security), as well as other potential numbers or reports. 
5 Further, use of the direct outputs can be made utilizing standard software tools, like Microsoft™ 
EXCEL™, to produce graphical representations of risk and value utilizing other models. Some 
examples of such are shown in Figures 2 through 5. These graphs were produced by inputting 
the values derived with the RAP methodology into Excel, with the desired values for the axes as 
shown, and running the graph routines having the shown data. As can be imagined, the reports, 

3o graphs, and other results are limited only by the requirements of the user. While particular 
elements, embodiments and applications of the present invention have been shown and 
described, it is understood that the invention is not limited thereto since modifications may be 

D made by those skilled in the art, particularly in light of the foregoing teaching. It is therefore 
contemplated by the appended claims to cover such modifications and incorporate those features 

|J5 that come within the spirit and scope of the invention. 
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What is claimed is: 

1 . An risk-adjusted pricing method for calculating the value of a derivative: 

a. determining the underlying security; 

b. determining the type of derivative; 

c. determining the risks associated with said derivative; 

d. determining the trading costs associated with said derivative; 

e. formulating the risk-adjusted pricing (RAP) equation for said derivative; 

f. solving said RAP equation for said derivative using numerical methods; and 

g. outputting a value for said derivative based on said solving of said RAP 
equation. 

2. An computer apparatus for calculating values of derivative securities comprising: 

a. an input unit taking an input of derivative characteristics and model parameters; 

b. a processing unit taking said input factors and computing a value for said 
derivative based on at least one of said inputs using a risk-adjusted pricing 
equation; and 

c. an output unit displaying said value of said derivative. 
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Abstract of the Disclosure 
A method and computer apparatus determines the value of a derivative by introducing 
risk premiums. The method includes determining the underlying security and derivative, 
determining the risks and trading costs associated with the derivative, formulating the RAP 
equation for the derivative, solving the RAP equation for the derivative and outputting a value of 
the derivative based on the solution of the RAP equation. The computer apparatus employs this 
method in a financial analysis application. The system includes an input unit for taking inputs of 
derivative characteristics and model parameters, a processing unit for computing a value of the 
derivative based on the input characteristics and parameters using a risk-adjusted pricing 
equation, and an output unit for displaying the value of the derivative. 
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# include <math.h> 
#include <stdlib.h> 
#include <iostream.h> 

#define MAX_TIME_STEPS 200 
#define MAX_PRICE_STEPS 200 
#define pi 3.1415962 
#define PUT 'P' 
#define CALL 'C 
#define MAX_GAMMA 128.0 

#define max(a / b) (((a) > (b) ) ? (a) : (b) ) 
#define min(a,b) (((a) < (b) ) ? (a) : (b) ) 

double 0[MAX_TIME_STEPS] [MAX„PRICE_STEPS] ; // array of the option price 
double Delta [MAX_TIME_STEPS] [MAX„PRICE_STEPS] ; // array of deltas 
double Gamma [MAX_TIME_STEPS] [MAX„PRICE_STEPS] 
double Theta [MAX_TIME_STEPS] [MAX_PRICE_STEPS] 
double Vega[MAX„TIME„STEPS] [MAX_PRICE_STEPS] ; 
double M[MAX_T1ME_STEPS] [MAX_PRICE_STEPS] ; // 
ments 
Int SO; 



// array of gammas 
// array of thetas 
// array of vegas 
array of the volatility adjust 



jtbid RAPoptionEuropean(// option's characteristics: 

double T, // T is time to maturity in years 
ffl double K, // K is the strike price in dollars 

W double S, // S is the current underlying price in dollars 

^ char Y, // Y is the option type: Put or Call 

ffl // model parameters: 



double sig, // 

double mu, // 

double R, // 

double C, // 

double // 



sig is the implied underlying volatility 
expected trend rate 

R is the variance risk premium coeff . 

c is the transaction costs/slippage coeff. 

v is the implied volatility risk premium coeff 



// DO 
double 
eter 
double 
double 
double 
double 
prices 
double 
double 
double 
double 
double 



double g, // gap risk coeff. 

double r, // risk free rate 

double rb, // borrowing fees rate 

double rd, // dividend yield rate 

// computational parameters: 

int nt, // number of time steps 

int ns, // number of price steps 

double Is, // lower bound on underlying price S 

double us // upper bound on underlying price S 

) { 

: initialize the computational variables 

ml = pow( 324.0 * R * C * C / pi, 1.0/3.0 ); // m is the risk param 



dt = T / nt; 
ds = ( us - Is ) / ns; 
s [MAX_PRICE„STEPS] ; 
s2 [MAX„PRICE_STEPS] ; 

t [MAX_TIME_STEPS] ; 
o [MAX_PRICE_STEPS] ; 
delta [MAX„PRICE_STEPS] 
gamma [MAX_PRICE_STEPS] 
theta [MAX_PRICE_STEPS ] 



// time step 

// the price step size 

// array of the underlying stock prices 

// array of squares of the underlying stock 



// 
// 
// 
// 
// 



array 
array 
array 
array 
array 



of the time steps 
of the option price 
of deltas 
of deltas 
of deltas 
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double vega[MAX_PRICE_STEPS] ; // 
double m[MAX_PRICE_STEPS] ; // 
double o_old[MAX_PRICE_STEPS] ; 



array of vegas 

array of the volatility adjustments 
// array of old option prices 



int it; //an auxiliary index variable 
int is; // an auxiliary index variable 
double ds2 = ds * ds; 
double a [MAX_PRICE_STEPS] 
double b[MAX_PRICE_STEPS] 
double c [MAX_PRICE_STEPS] 
double d[MAX_PRICE_STEPS] 



double auxl, aux2, aux3 , aux4, aux5, aux6; 



// array of matrix coefficients 

// array of matrix coefficients 

// array of matrix coefficients 

// array of matrix coefficients 



// 
if 
if 
if 
if 
if 
if 
if 

.v=e 
!i? v 

i^pro 



DO: Make sure that the coefficients are within their meaningfull range 



( nt > MAX_TIME_STEPS 
( ns > MAX_PRICE_STEPS 
( sig < 0.0 ) sig 
( R < 0.0 ) R = -R; 
( g < 0.0 ) g = -g; 
( rb < 0 .0 ) rb = -rb; 
( Is < 0.0 ) Is = 0.0; 



( us < 0.0 



us = -us; 



nt < 1 ) nt = MAX_TIME_STEPS ; 
ns < 1 ) ns = MAX_PRICE_STEPS ; 
sig; // volatility cannot be negative 

// variance risk premium cannot be negative 
// gap risk premium cannot be negative 
// borrowing fees cannot be nagative 

// underlying price range cannot be negati 

// underlying price range cannot be negati 



=0.0; // set the volatility risk premium to zero in this demonstration 
cedure 



M // DO: Set the borrowing rate and dividend rate for puts to zero 
ffl if ( Y == PUT ) { 
M 1 rb = 0.0; 

L } 

jE // DO: Initialize the arrays of the time 
5 for( it = 0; it <= nt; it++ ){ 

U t[it] = it * dt; // the maturity time T is t[nt], the current time 0 is 
lfc[0] 

i } 

// DO: Initialize the arrays of the price 
is = 0; 

s[0] = S - is * ds; 
while ( s [0] >= Is ) { 

s[0] = s[0] - ds; 

is++; 

} 

SO = is - 1; 
s[0] = s[0] + ds; 
for( is = 1; is <= ns; is++ ){ 
s[is] = s[is-l] + ds; 

} 

for( is = 0; is <= ns; is++ ){ 

s2[is] = s[is] * s[is]; // array of the price squares 

} 

// DO: initialize the finite difference scheme by the derivative's values 
/ / at maturity time by the payoff function 

if ( Y == CALL ) { 

for( is = 0; is <= ns; is++ ){ 

o[is] = max( 0.0, s[is] - K ); // call option payodeltasff 
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if ( s[is] < K ) { 

delta [is] =0.0; // delta zero below the strike 

gamma [is] =0.0; // gamma zero 

} else if ( s[is] == K ) { 
delta [is] =0.5; //at the strike delta 0.5 for numerical stability 

gamma [is] = MAX_GAMMA; // gamma is the maximum allowed by the market 

// which is inverse of the minimum tick size 

} else if ( s[is] > K ) { 
delta [is] =1.0; // delta one above the strike 
gamma [is] =0.0; // gamma zero 

} 

theta[is] =0.0; //no time decay left 
vega[is] =0.0; //no volatility sensitivity 

} 

} else if ( Y == PUT ) { 

for( is = 0; is <= ns; is++ ){ 

o[is] = max( 0.0, K - s[is] ); // put option payoff 

if ( s[is] < K ) { 
delta [is] = -1.0; // delta zero below the strike 
gamma [is] =0.0; // gamma zero 

} else if ( s[is] == K ) { , 
delta [is] =0.5; //at the strike delta 0.5 for numerical stability 
gamma [is] = MAX_GAMMA ; // gamma is the maximum allowed by the market 

// which is inverse of the minimum tick size 

} else if ( s[is] > K ) { 
delta [is] =0.0; // delta one above the strike 
gamma [is] =0.0; // gamma zero 

theta[is] =0.0; //no time decay left 
vega[is] =0.0; //no volatility sensitivity 

} 

} 

for( is = 0; is <= ns; is++ ){ 

if( gamma [is] >0.0){ r • •, , n /-o n n 

m[is] = max( 0.0, 1.0 + g - min{ ml * pow( gamma [is], 1.0/3.0 ), 

2.0 * R * delta [is] * delta [is] / gamma [is] ) ); 

} else { 

m[is] = 1.0 + g; 

} 

} 

for( is = 0; is <= ns; is++ ){ 
0[nt] [is] = o [is] ; 

Delta [nt] [is] = delta [is]; 
Gamma [nt] [is] = gamma [is]; 
Theta[nt] [is] = theta[is]; 
Vega [nt] [is] =vega[is]; 

} 

// DO: remember the last computed values 
for( is = 0; is <= ns; is++ ){ 
o_old[is] = o[is] ; 

} 

// DO: Setting: starting time is maturity time 
it = nt; 

//DO: set the auxiliary variables for faster computational speed 
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auxl = 0.25 * sig * sig * dt / ( ds * ds ) ; 
aux2 = 0.25 * ( r - rb - rd ) * dt / ds; 
aux3 = 0.5 * r * dt; 

// DO: in the loop compute all the optionis prices 
while ( it >= 0 ){ // INV: new time >= current time 

it--; // Meaning: Time := Time - Time Step 

//DO: compute the solution on the new time level by solving 

// the RAP equation with the semi-implicit finite difference 

// Crank-Nicolson scheme method 

//DO: Initialize the triangular matrix coefficients and the right hand 

side 

for( is = 1; is < ns; is++ ){ 
// a,b,c are on the diagonal 
aux4 = auxl * s2[is] * m[is] ; 
aux5 = aux2 * s [is] ; 

a [is] = aux4 - aux5; 

b[is] = -1.0 - 4 * aux4 - aux3 ; 

4f c[is] = aux4 + aux5; 

* y // d is the right hand side 

J d[is] = o_old[is-l] * ( - aux4 - aux5 ) + o_old[is] * ( -1 + 4 * aux4 
"f^ : a ux 3 ) 

W + o_old[is+l] * ( - aux4 + aux5 ) + v * vega[is]; 



// DO: Initialize the boundary conditions: 
if ( Y == CALL ) { 

// Call boundary conditions 

b[0] = 1.0; 

c[0] = 0.0; 

d[0] = 0.0; 

a[ns] = 0.0; 

b[ns] = 1.0; r . n xx 

d[ns] = ( s[ns] - K ) * exp ( ( r - rd - rb ) * ( T - t[it] )); 
} else if ( Y == PUT ) { 

// Put boundary conditions 
b[0] = 1.0; 
c[0] = 0.0; 

d[0] = ( K - s[ns] ) * exp( - r * ( T - t[it] )); 

a[ns] =0.0; 
b[ns] = 1.0; 
d[ns] = 0.0; 

} 

// DO: Solve the matrix for the new options values 
// DO: run the forward elimination first 
for( is = 1; is < ns; is++ ){ 

aux6 = - a [is] / b[is-l]; 

b[is] += c[is-l] * aux6; 

d[is] += d[is-l] * aux6; 

} 

// DO: run backwards elimination to get the options values 
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o [ns] = d[ns] ; 

for( is = ns-1; is > 0; is-- ){ 

o[is] = ( d[is] - c[is] * o[is+l] ) / b[is]; 

} 

o[0] = d[0] ; 

// DO: compute the sensitivities 
for( is = 1; is < ns; is++ ){ 

o[is] = max( 0.0, o[is] ); // option value should not be negative 

delta[is] = ( o[is+l] - o[is-l] ) / { ds + ds ) ; 

gamma [is] - min( MAX_GAMMA, 

max( 0.0, ( o[is+l] - 2 * o[is] + o[is-l] ) / ds2 ) ); 

theta[is] = ( o_old[is] - o[is] ) / dt; 

vega[is] =0.0; //in this software vega need not be computed 

// because we assume v zero 

if ( gamma [is] > 0.0 ) { 
m[is] = max( 0.0, 1.0 + g - min( ml * pow( gamma [is], 1.0/3.0 ), 

2.0 * R * delta[is] * delta[is] / gamma[is] ) ); 

} else { 
m[is] = 1.0 + g; 
} 

// DO: save the previous step for the options 
o_old [is] = o [is] ; 

} 

//DO: Save the computed values into the arrays 
for( is = 0; is <= ns; is++ ){ 

0[it] [is] = o[is] ; 

Delta[it] [is] = delta[is]; 

Gamma[it] [is] = gamma[is]; 

Theta[it] [is] = theta[is]; 

Vega [it] [is] = vega [is]; 

M[it] [is] = m[is] ; 

} 

} 

// INV: the values have been computed up to the current time 

// Return of the computed derivative's values and sensitivities 

//is done throughthe globally defined arrays 0, Delta, Gairnna, Theta, Vega 
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Options 



No mystery 
behind 
the smile 

Milan Kratka presents a risk-adjusted pricing model that 
matches the market prices of options more closely 




lassical option theory derives option values from the 
so-called arbitrage replication process: investors who 
buy options and hedge them continuously and with- 
out cost are expected to generate a risk-free return 
equal to that from Treasury bonds. These idealistic 
assumptions imply that the option model is not con- 
sent with real markets. In real life, every investor rehedges at discrete 
tim| intervals, pays trading costs, is exposed to potential price gap and 
volatility risks, and experiences positive portfolio variance. 

fjhe pricing model we present here helps close this gap between theory 
afilf practice. Risk-adjusted pricing incorporates all the real-life factors into 
a -simple but comprehensive model that is consistent with the market. It 
esfences the old "no-arbitrage" paradigm with the more realistic notion 
oFSalue-at-risk. We will try to explain the main ideas behind this break- 
through methodology. 

Standard assumptions 

Suppose that an investor prices an option 0 on a stock S that pays zero 
dividends and follows the typical, constant-volatility lognormal process: 

dS = jtSdt + cSdl 

where \i denotes the drift (or the trend) and 0 denotes the volatility. Let 
us assume we buy N options and short-sell 6N stocks as a hedge. We plan 
to rehedge this portfolio by trading in stock. Let us denote the rehedging 
time interval by A t . This should reflect the trading strategy chosen by the 
investor and, in general, will vary with the time to expiry, underlying price 
and other factors. Since no portfolio can be hedged continuously or with- 
out cost, the investor needs to specify a rehedging rule. For example, we 
might simply want to hedge in the most advantageous way. So we ask 
our quants to figure out how to do this. They may come up with strate- 
gies such as: 

□ rehedge on sunny days only; 

□ rehedge every hour; 

□ rehedge with every 1% change in the stock price; 

□ rehedge with every 0.05 change in delta; or 

□ leave rehedging strategy to the best trader. 

What would be the optimal strategy? Rehedging too frequently accu- 
mulates high trading costs. Hedging too infrequently results in high port- 
folio variance. The risk-adjusted pricing methodology offers a natural 
solution to the problem of specifying an optimal rehedging strategy. 

Risk-adjusted pricing 

We have now set up a portfolio that will be rehedged in a time interval A t . 
Over this period, we expect it to generate a return equal to that from risk- 
free securities, but it should also pay for the trading costs and compensate 
us for having positive portfolio variance and for taking volatility exposure 
risk. Though the transactions will occur at discrete times, we can locally 



estimate the rate of trading costs by spreading them over the time between 
the rehedging points, as is the customary practice with the dividend rate. 
We will estimate the rate of the portfolio variance similarly. 

Transaction fees, trading strategies and risk aversion differ from one 
investor to another. An option may seem to be undervalued for one in- 
vestor and overpriced for another. The market price of a liquid securi- 
ty reflects the prevailing expectations among all active investors. Who 
would not like to have a pricing model which is both in line with the 
market and reflects investor-specific conditions and expectations? The 
latter is not quite possible without the former. Risk-adjusted pricing al- 
lows us to model both. 

Let us denote the risk-free rate by r, the stock borrowing rate by r B , the 
expected rate of trading costs by r TC? the expected rate of premium for tak- 
ing positive variance risk by r VAR and the expected rate of premium for 
volatility exposure risk by r a . The average change E[dP] of the portfolio 
value in an infinitesimal time increment dt should then satisfy: 

E[dP] = E[dO]-E[5dS] 

= (9 t 0 + |a 2 S 2 9 S sO)dt + E[(3 S 0 - 8)dS] 

= r(0 - 5S)dt + r B SSdt + (r VA R + r TC + r c )dt 

Rehedging to the market delta 3 S 0 would make the expected value 
E[{d s 0 - 8)dS] equal zero, and the portfolio would have zero variance over 
the time dt. On the other hand, rehedging consistently to a different delta 
5 would make the expectation term equal to (3 S 0 - 8)p,Sdt ) and it would 
increase the expected portfolio variance by 0 S O - 5) 2 o 2 S 2 dt. For investors 
with strong trend expectations, it is optimal to trade to a trend-adjusted 
delta that differs from the market delta. 

In general, it is difficult to estimate the future trend term from the 
market. It is a subjective issue. If the trend was known and had significant 
value, the underlying price would have already been adjusted according- 
ly. There is much uncertainty about trends. Maybe that is why they are 
often referred to by the more moderate term "drifts". Each time we buy a 
stock, somebody sells it. Similarly, some investors might underhedge, while 
others might overhedge. For the market modelling we could simply as- 
sume that most traders hedge to the market delta. For investors with a 
strong opinion about the trend, or even for a market reflecting these opin- 
ions, rehedging to a trend-adjusted delta should be more advantageous, 
since it should increase expected return. 

Portfolio variance risk premium 

How much, in premium, should investors be compensated for positive 
portfolio variance risk? It is safe to say that the higher the variance, the 
higher the premium. As the variance of a normal process increases linearly 
in time, it may seem natural to expect the premium to follow the portfo- 
lio variance linearly. This is similar to portfolio performance rating using 
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Sharpe Ratios. If we denote the variance risk coefficient by R, the rate of 
the portfolio variance risk premium is* 

R - var(AP) 

You can also consider R as a risk aversion coefficient. The variance of 
the portfolio between two rehedging points can be estimated as: 

"(AP-E[AP]f 



var(AP) = E 



(3 s O(t + A t )-5(t)) 4 
(AS) 2 



r 2 (AS) 4 + (a s o(t)-8(t)) 2 (AS) 2 



= 3r 2 a 4 S 4 A t 2 + (3 S 0 - 5)Vs 2 A t 

where AS denotes the change in the underlying price between rehedging 
times t and t+A t , and T = d ss 0 denotes market gamma, as is usual The 
rate of the portfolio variance risk premium is then: 

% r VA R = 3Rr 2 a 4 S 4 A t + R (3 S 0 - 5) Vs 2 

J^otably, the longer the time between rehedging A t , the higher the rate 
olWie variance nsk premium r VAR . 

fiior a risk-averse investor with a strong opinion about the trend, the op- 
XmM delta to rehedge to is the one that maximises the risk-adjusted, trend- 
related rate: 



(3 S 0 - 8)(m- - r + r B ) S - R(d s O - 8) Vs 2 



j^hat delta is. 



8 = 3 s O- 



2Ra 2 S 



ssSThis trend-adjusted delta makes the risk-adjusted, trend-related rate 
(yfr r + r B ) 2 /4Ra 2 A moderate up-trend would cause option holders to 
b'fjiedged with a lower delta than the market delta. Trends seem to have 
amatively high impact on options that are not at-the-money. Whatev- 
er the direction and strength of future trends, their presence and mar- 
ket anticipation will help to enlarge the effect of higher implied volatility 
at the tails. 

Trading costs estimator 

An investor's trading costs typically depend on the number of contracts 
traded, type of order, size of bid-ask spread, bid and ask volumes, ex- 
change fees and a transaction fees structure negotiated with a broker. Let 
us denote an aggregate trading cost per contract, when the size of a trade 
is n contracts, by C{n). It is typically a non-increasing function of a trade 
size, C(n) < 0. The larger the trading size, the higher the total commissions 
but the lower the total cost per traded share. It does not need to be posi- 
tive, as heavy buying may move the market above the average cost of the 
trade, causing an immediate paper gain. For market modelling, we may 
assume that the prevailing aggregate trading cost per contract is constant 
C(n) = C. Investors should specify their own trading costs structure. It will 
affect investors' expected option theoretical values and their choice of re- 
hedging strategy. Trading costs alone would have a relatively small price 
effect but combined with discrete hedging, positive variance risk and volatil- 
ity exposure premiums can have a significant impact. 
The expected rate of trading costs: 



lTC = 



C(N|A5|)|A5| 



reflects that we will trade N I A5 1 shares and spread out the cost over a time 
interval A t . The average amount of shares traded per option in each delta 
rehedge adjustment can be estimated using Ito's lemma: 



E[ | A5 1 ] = E a t SA t + a S 5AS + |d S sS(ASf 
= aS|d s 5|E[|dz|] 



If the time between rehedging was long enough, the expected value 
could be calculated more precisely, and would marginally depend on the 
trend \i. The rate of trading costs can finally be expressed as. 

rTc^^cslrlV^^oSlrl-^ 

Notably, the longer the time between rehedging A t , the lower the trad- 
ing costs rate r TC . 

The higher the volatility sensitivity of the portfolio (w T hich equals the 
option's vega), the higher the volatility risk we are exposed to. How should 
we model the volatility risk exposure? It will depend on our anticipation 
of volatility change, which is a subject for a separate project, a volatility 
model We can model the volatility risk premium by r 0 = vl d a 0j . sa\ :e- 
flecting our portfolio sensitivity to any change in underlying volauhn The 
higher the vega, the higher the premium. The specific form of the volatil- 
ity risk premium coefficient should come from our volatility model. It should 
reflect the probability of the volatility change as well as the size of its prob- 
able change. 

Optimal hedging strategies 

If we believe in the efficiency of financial markets, we might suppose that 
market option prices reflect the market value of risk, average prevailing trad- 
ing costs and competitive rehedging strategies. There is never a guarantee 
that the free market will value options in any particular way. But we can as- 
sume it does, build a pricing model and use it to verify our presumption 

Given investor-specific trading costs C(n) and an investor-specific vari- 
ance risk coefficient R, the optimal rehedging strategy is identified by min- 
imising the sum of r VAR and r TC , the variance premium and the trading cost 
terms: 

r V AR + r T c = 3Ra 4 S 4 r 2 A t + C No S| T | ^j^cS\ T | 

where the rehedging time A t varies in a reasonable range. This may span 
a fraction of an hour to a fraction of a year. The gamma, T, should be equal 
to the market gamma d s 0. The simple task of finding the minimum can be 
done approximately or even precisely. There is no need to write elaborate 
papers about what an optimal hedging strategy in the presence of trading 
costs ought to be. If the trading costs per share, C, do not depend on trade 
size, the optimal rehedging time would be: 

A t = - 



a 2 S 2 r 2/3 



where we have denoted the constant: 



k = (C/3RV^)' 



1/3 



This translates to a very simple strategy of rehedging with every stock 
change of: 

as = +<jsVa7 = ±kr~ 1/3 

This is probably close to a rehedging rule you have already used. If not, 
feel free to upgrade now! 

In line with our intuition, we shall trade more often with higher risk 
aversion R, lower trading costs C or higher gamma F. In the limiting case 
of zero transaction costs, we would trade continuously, as assumed by the 
Black-Scholes model. 

The assumption of competitive hedging in the market will then yield 
the sum of the variance premium and the trading costs terms to be: 

rvAR = r TC = ^ma 2 S 2 rr 1/3 
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1. Black-Scholes implied volatility 




where we have denoted the constant m = (324RC 2 /jt) 3/3 . This was the last 
-:.@ece needed to finalise a surprisingly simple risk-adjusted pricing equation: 

dtO + \ <^S 2 {l - mr 1/3 )r = r(0 - SS) + r B 6S 

It can be made more complete by including the drift term -(u- r 4- r B ) 2 /4R<£ 
and the volatility risk premium term that in our simplified form was v I d ff 0 1 . 
It can be solved numerically backwards in time by finite difference scheme 
methods. The three-dimensional sample plots above (figures 1-4) have been 
produced by choosing an underlying volatility of 30%, a risk parameter of 
0.9 and an interest rate of 6%. For the purposes of this paper, we have left 
out the technical details of the model implementation. 1 We will just note 
that the rate of the portfolio variance should not be greater than the vari- 
ance 5 2 c 2 S 2 of an unhedged option position, which caps the variance risk 
premium that would have otherwise spiked up at the strike price at the op- 
tion's maturity. 

Volatility smile 

There are just two market parameters in the simple risk-adjusted option 
pricing equation - the anticipated volatility of the underlying price, o, and 
the so-called risk parameter m - although we would consider the volatil- 
ity risk premium coefficient v, trend parameter \i and risk aversion coeffi- 
cient R for a more detailed option model. The volatility risk premium 
coefficient measures our exposure to the risk that the market could change 
its mind about the anticipated volatility. Its effect on option pricing should 
diminish with time to maturity, and thus it would not explain sharpening 
of the so-called volatility smile. The risk parameter is the one responsible 
for its sharpening, as its effect on option pricing magnifies with gamma, 
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2. Risk-adjusted call price minus Black-Scholes 
call price 




4. Risk-adjusted gamma minus Black-Scholes 
gamma 




as if the options with higher gamma were priced with lower volatility. The 
trendiness of the market ji helps to sharpen the smile at the tails. 

In any case, it is quite nonsensical to assume a different underlying 
volatility for different options. There can be only one anticipated under- 
lying volatility, no matter which options we look at. Even if we could in- 
corporate fat tails, price gaps, complex supply-demand relations or even 
autocorrelation, there would still be only one underlying process with only 
one volatility. The notion of implied volatilities points out how far away 
the market pricing is from an idealised risk-neutral Black-Scholes world. 

So, why do we observe volatility smiles? Because there are costs and 
risks associated with options trading, as we have seen and become capa- 
ble of modelling. There is no mystery behind the smile any more. The risk- 
adjusted pricing methodology may not replicate the market pricing 
perfectly, but it does explain quite a significant part of its stable dynam- 
ics. It provides an edge in understanding option valuation and optimal 
hedging strategies. It shows how to incorporate various market risks quan- 
titatively into an option pricing equation. It will even save resources that 
would otherwise be spent on artificial implied volatility surfaces model- 
ling. Finally, it is simple and easy to put into production. 

We have chosen to model market risks in a simplified way to demon- 
strate the natural power of the risk-adjusted pricing methodology. It is just 
a first Step in pointing out the new dimensions possible in the financial 
modelling of a not quite risk-neutral world. ■ 
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